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Abstract 

The stability of the quenched incommensurate phase in two dimensions 

against the creation of overhangs and finite loops (OH/FL) in the replica 

space is investigated for a model of domain walls with N colors. Introducing 

a chemical potential e for OH/FL, the probability for the formation of these 

objects is studied for e — > 0. In the pure limit this probability vanishes with 

y/e, whereas the fluctuations of this probability are long-range correlated in 

the quenched system. This indicates an instability related to the symmetry 

in replica space. It is accompanied by the creation of a massless boson. The 

latter leads to a power law decay with exponent oc 1 /N for the product of the 

correlation functions along the domain walls. 
PACS numbers: 64.70.Rh, 64.60.Cn, 64.60.Fr 
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The model for the commensurate-incommensurate transition (CIT) in d = 1 + 1 dimen- 
sions is the prototype for a class of systems characterized by directed domain walls, directed 
random walks, directed polymers or flux lines. The common feature is that of interacting 
random walks where the walkers choose randomly steps forward, left or right but not back- 
ward. The interaction is due to the condition that walks are not allowed to cross. In the 
presence of quenched disorder these systems should exhibit a new phase which is related to 
a glass or frozen phase. The latter can be regarded as freezing of the domain walls (directed 
random walks etc.) in a random potential. From a conceptional point of view the CIT 
model has attracted attention because it is soluble without disorder Q, and the effect of 
disorder can be regarded as a perturbation. The disordered case has been studied by various 
methods like the Bethe ansatz @], scaling arguments ||, variational approaches [§J and with 
perturbation theory ||. Although these works agree on the fact that the critical exponent 
of the density changes from 1/2 in the pure system to 1 in the disordered system they came 
to different conclusions concerning the decay of the correlation function in the quenched 
system. In a recent article by Tsvelik || it was argued that quenched disorder is irrelevant 
with respect to the asymptotic decay of the density density function. This is surprising 
since the correlation function should be more sensitive to additional fluctuations than the 
density. Therefore, at least near the CIT one should expect that disorder changes the qual- 
itative properties of the model. For the corresponding model with disorder correlated along 
the direction of the domain walls was found that the long-range correlation is destroyed 
0. In the following we will present a calculation which shows that there is indeed a new 
feature due to disorder which affects the correlations. It is characterized by the formation of 
overhangs and finite loops (OH/FL). A crucial point in the various approaches to the CIT 
in the presence of disorder is the replica trick. It has been argued that there is no replica 
symmetry breaking (however, see ||). In contrast to this we will start from a broken replica 
symmetry given by an external field (chemical potential) which creates OH/FL. 

The grand canonical statistics of domain walls in d = 1 + 1 dimensions can be described 
by a fermion Lagrangian density in imaginary time representation || 
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L = c + d T c — d x c + d x c — v(x, r)c + c + /ic + c (1) 
and the partition function 

//■oo r- 
P[c + ,c]exp(-y o dr J dxL). (2) 

The imaginary time r is along the direction of the domain walls and x is perpendicular to r. 
Disorder is introduced by the random potential v(x, r) which is Gaussian distributed with 
zero mean and (v(x,t)v(x',t')) = a 2 5(x — x')5{t — r'). 

The partition function is invariant under a time-reversal transformation. Only the time 
differential operator is changed as d T — > = —d T by this transformation. The integration 
of the non-interacting fermions in Z gives the space-time determinant det(<9 T + d 2 + \l — v). 
Then the time-reversal invariance of Z is obvious because the transposition of the matrix 
leaves the determinant invariant. 

Averaging of the free energy or of an expectation value can be performed either by using 
a supersymmetric generalization of the fermion Lagrangian || or by using the replica trick 
InZ = lim n _ > o(^ n — l)/ n - Since both representations are equivalent if we choose the correct 
structure, we will use the replica trick here for convenience. Introducing an even number of 
replicas, say 2n, we take one sector of n replicas with d T and the other sector of n replicas 
with &l ' . This can be written in a spinor representation for the Lagrangian density as 

d+) \d T T+ dl-v + ^ ) \c 

where c, d are n-component fermions, and the matrix in (|3|) is diagonal w.r.t. to the n 
replica components. Now we introduce a chemical potential e which allows the formation 
of OH/FL by combining the two different replica sectors. This means that the chemical 
potential creates a particle-hole pair c + d which can be annihilated at another space-time 
point by d + c. Including this new chemical potential in the Lagrangian density (§) we replace 
the zeros in the diagonal elements by ie. The imaginary unit % guarantees a positive weight for 
the OH/FL in the partition function because OH/FL contain always a sequence of reversed 
Grassmann variables which contribute a minus sign. The density of OH/FL can be measured 
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by varying locally the chemical potential: p e (x,r) oc (c + (x, r)d(x, r) + d + (x,r)c(x,T)). In 
the pure limit v = this density can be evaluated by a simple calculation as p e ~ const. y^e; 
i.e., it vanishes continuously for vanishing chemical potential e. Thus, the pure system is 
stable against the formation of OH/FL. Another interesting quantity is related to the fact 
that OH/FL are given by a combination of local edges (s. Fig.l) which are separated by a 
distance in x and r. The correlation of an edge at (x, r) and another one at (0, 0) is 

(c + {x, r)d{x, r)d + (0, 0)c(0, 0)). (4) 

Because this is a correlation function of non-interacting fermions it will factorize into 
(c + (x,T)d{x,T))(d+(0,0)c{0,0)) - (c + (x,T)c{0,0))(d + {0,0)d(x,T)). The first term measures 
the density of OH/FL whereas the second term is the product of correlations along the 
domain walls from (0, 0) to (x, r). 

Our calculation will be based on the generalization of the domain wall model to one which 
has walls with N different colors. Without quenched disorder the new model separates 
into N independent models. However, the disorder potential is chosen as v a p{x,r) with 
a, (5 — 1,...,N such that the domain walls may change statistically the color from a to 
p. This mechanism reduces the hard-core of the domain walls (i.e., Pauli's principle of the 
fermions) because the domain walls can avoid each other by changing the color. The complex 
matrix elements v a p are statistically independent except for the symmetry v a p = Vp a . After 
averaging we obtain an effective interaction for the fermions (c = (d, c)) in the replica model 

E c+A^c+jC^. (5) 

i,j=l a,(3=l 

The limit iV — > oo can be solved by a saddle approximation. This becomes obvious if one 
decouples the effective fermion-fermion interaction of the quenched model by a new random 
matrix field Q which includes fluctuations of both chemical potentials \i and e: we replace 
a 2 N- 1 c+ a Q i/3 c+ 3 c jiQ! = -^N^c^Cj^c+pCi^ by 2iQ ji J2 a cf, a Cj,a- Q is a Hermitean 2n x 2n 
matrix field. Going back to the replicated partition function Z 2n we integrate out the fermion 
field c ijQ . This leads to an effective action which depends only on the decoupling field as 
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r°° r 2 
S e ff = N J Q dr J dx—Tr 2n (Q(x,T) 2 ) 



it + 2iQ 2 \ (x, r) d T + d 2 + \x + 2iQ u (x, r) \ 

iVlndet | (6) 
+ d% + u. + 2iQ 22 {x,r) ie + 2iQ 12 (x,r) J 

Tr 2n denotes the trace w.r.t. the 2n replica and det the determinant w.r.t. the 2n replica and 

x, t. This is the 1+1-dimensional replica version of the 2 + 1-dimensional 'supersymmetric' 

action found for the flux lines in a random potential . 

S e ff depends on the number of colors only through the prefactor N. Therefore, the 

effective field theory for large N enables us to do the functional integration in saddle point 

approximation. This approximation can be interpreted as the replacement of the chemical 

potentials \x and e in the free fermion propagator by a self-energy term which is determined 

by the saddle point equation SS e ff = 0. There are two contributions if we assume that 

Q provides only additive corrections to the two chemical potentials. One is a shift of the 

chemical potential /x in the limit e — > 

1 2 

/i s = 2iQ jj = -a 2 J dkdu (J = !>•••> 2 ") (7) 

with D(e') = e' 2 + uj 2 + (— k 2 + fi + fi s ) 2 , and the other is a shift of the chemical potential 
for OH/FL 

e' = 2Q jj+n = 2Q j+nj = eV J dkdco/D(e r ) (j = 1, .., n). (8) 

e' = is always a solution of (|8|). This is a replica-symmetric (RS) solution. A replica- 
symmetry breaking (RSB) solution e' ^ can be found if 

J dkdu/D(0) > 1/a 2 . (9) 

Since the denominator D(e') is increasing with increasing e' 2 there is an e' which satisfies 

/" dkdu/D{e') = 1/a 2 . (10) 

This can be rewritten as 
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l+u 2 + {-k 2 + n'/e') 2 ' 



In this case we obtain from (|7J) a renormalized chemical potential 




(12) 



-D(e') is increasing with decreasing p! < 0. This implies that there is a p c < (which 
depends on a) such that J dkdu/D(0) < 1/cr 2 for //' < /i c . Consequently, there is only a RS 
solution e' = for p! < p c . This behavior describes a transition from a RS solution to a 
RSB solution if we change p or the disorder a. 

The evaluation of the density of domain walls in saddle point approximation requires a 
cut-off, because only a finite number of domain walls per unit length is reasonable. With 
the cut-off k 2 = 1 we obtain after integration over uj 



p! < and e' — > implies p — > 0. The density vanishes linearly. For instance, the density of 
the RS solution is p oc o 2 + fi. Thus the density of the RS as well as the RSB solution are 
in agreement with the linear behavior of the Bethe ansatz calculation of Kardar 0. 

Finally we evaluate the Gaussian fluctuations around the saddle point solution which are 
the 1/N corrections. The stability matrix can be taken from Ref. 0. The fluctuations of 
Qn, Q22 are massive with the eigenvalues A± 



i.e., the eigenvalues of the RSB solution are always positive (stable). On the other hand, for 
e' = the eigenvalue A„ is 




(13) 




(14) 



For e'^0A_ is 




(15) 




(16) 



6 



Since inequality (||) holds if there are two saddle point solutions, A_ of the RS solution 



(|16|) is negative (unstable) when the RSB solution exists. The other two eigenvalues of 
the RSB solution are A 3 = 2(A + — A_) > (related to Qu — Q21) an d the massless one 
A4 = (related to Q12 + Q2i)- The massless fluctuations are a consequence of a global 
symmetry of the model 0. They can be expressed as a non-linear sigma model for Q with 
the constraint Q 2 = 1, Tr2 n Q = and Q02Q = — The constraint can be satisfied if we 
parametrize Q by an n-component field ip as Q = 03 (cos 2 p — sin 2 ip) + 2o\ cos p sin p. This 
parametrization neglects a unitary transformation inside of each replica sector which is not 
expected to affect the properties of the ^-dependent part of the model. With the field p the 
action of the non-linear sigma model reads 

S eff = bN J dxdrp(d 2 x + d 2 )p (17) 

with a positive constant b, depending on 7 but not on iV or e, which are given in Ref. 0. 
From this bilinear action we can evaluate the correlation function of the density of the edges: 

(Qi2(z,r)Q 21 (0,0)) ~e' 2 const. (x 2 + r 2 y 1/2TbN (18) 

for large distances. Thus there is no replica symmetry breaking because the correlation of 
the order parameter decays. According to our previous remark this implies also a power 
law decay for the product of correlation functions along the domain walls. In contrast, a 
single averaged correlation function decays exponentially according to the finite length scale 
created by e'~ l . On a pure surface the correlation function decays like (x 2 + r 2 ) -1 / 2 . 

In conclusion, we have found that the statistics of domain walls in the CIT model with 
quenched disorder is unstable against the formation of OH/FL in the replica space. This 
result was based on a saddle point calculation for an effective self-energy in the replica 
model with 2n replica and iV colors. The properties of the model are similar to those 
of the Gross Neveu model |Hj and other 1+1- dimensional fermion models with attractive 



interaction [pT| . The creation of OH/FL in the quenched CIT model is plausible because 



the random potential favors OH/FL which can easily circumvent an unfavorable potential 



or freeze into the local potential wells. Since there are no OH/FL by definition, the replica 
model creates them spontaneously by combining wall elements coming from different replica 
sectors. However, the effect of the OH/FL is only marginal in our calculation because 
they appear with vanishing density. Nevertheless, the instability the instability related to 
the OH/FL means a drastic change in the quenched state compared with the pure CIT 
problem. 

A similar instability was also discussed in the context of the Gross Neveu model by 
Witten [12 1 who evaluated the correlation function of the massless bosons. Those bosons 
also lead to a power law decay with an exponent proportional to 1/N. 

It should be emphasized that the formation of OH/FL is a special type of instability 
regarding the freezing into local minima of the random potential. We cannot exclude other 
instabilities which may be even more favorable. Therefore, the purpose of this investigation 
was only to provide an example for an instability in the quenched CIT model. It indicates 
that the model has probably a rich structure which deserves further investigations based on 
effects related to replica symmetry. 
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FIGURES 

FIG. 1. A finite loop (a) and an overhang (b) created for the domain wall segments from 
different sectors of the replica space j and k. 
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